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Abstract 

We investigate the spacetime of anisotropic stars admitting conformal motion. 
The Einstein field equations are solved using different ansatz of the surface tension. 
In this investigation, we study two cases in details with the anisotropy as: [1] 
Pt = npr [2] pt — Pr = g^Cf? + C2) where, n, ci and C2 are arbitrary constants. The 
solutions yield expressions of the physical quantities like pressure gradients and the 
mass. 
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1 Introduction 



Most of the stars in the galaxies are the main sequence stars which evolve by burning 
lighter elements into heavier nuclei. Stars massive than ~ lOM© explode into supernova 
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leaving the core behind which then collapses to form a compact object. The cores are 
supported by the degenerate pressure of its constituent particles and possess the densities 
of the relativistic scales i.e. Rg = 2GM/rc^ ~ 1 (which is the Schwarzschild radius of 
the star) where M is the mass and r is radius of the compact object. At these densities, 
relativistic effects dominate and the physical quantities like gradients of pressure and mass 
are determined by the Tolman-Oppenheimer-Volkoff (TOV) equations for an isotropic 
and homogeneous compact star. In these stars, the matter is found in stable ground 
state where the quarks are confined inside the hadrons. It is suggested that the quarks 
if de-confined into individual u, d and s quarks can also yield a stable ground state of 
matter, which is called 'strange matter' [1]. Stars composed of mostly strange matter 
are therefore termed as 'strange stars' [2] (see [3, 1] for review on this topic). Generally, 
superdense stars with mass to size ratio exceeding 0.3 are expected to be composed of 
strange matter [5]. The prime motivation for the existence of strange stars was to explain 
the exotic phenomena of gamma ray bursts and soft gamma ray repeaters [6, 7]. Now 
with the observations of the Rossi X-ray Timing Explorer, it is convincingly shown that 
astrophysical source SAX J1808. 4-3658 is more likely a strange star [8]. The transition 
from the normal hadronic to the strange matter occurs at sufficiently high densities or 
corresponding low temperatures as T oc oc p~^. These conditions can mostly be found 
inside the cores of fast rotating pulsars, P-stars (composed of u and d quarks and are in /? 
equilibrium with electrons) or magnetic field powered magnetars. The density distribution 
inside these stars need not be isotropic and homogeneous (as proposed in the TOV model) 
if strange matter truly exists, then stars composed of entirely strange matter can also 
be found. Recently several authors have studied compact stars with anisotropic matter 
distribution [9, 10, 11, 12, 13]. We provide star models admitting conformal motion with 
different anisotropy. 

To search the natural relation between geometry and matter through the Einstein's Equa- 
tions, it is useful to use inheritance symmetry. The well known inheritance symmetry is 
the symmetry under conformal killing vectors (CKV) i.e. 

LiQik = ipgik, (1) 

where L is the Lie derivative operator, ^ is the four vector along which the derivative is 
taken and ip is the conformal killing vector. Note that ii ip = then Eq. (1) gives the 
Killing vector, if ip = constant it gives homothetic vector and if ip = ip{'x, t) then it yields 
conformal vectors. Thus CKV provides a deeper insight into the spacetime geometry. 
Moreover, if the conformal factor ^/; = 0, it implies that the underlying spacetime is con- 
formally fiat which further implies that the Weyl tensor also vanishes. These conformally 
fiat spacetimes represent gravitational fields without the source of matter producing these 
fields. 

The plan of the paper is as follows: In the second section, we model our gravitational 
system and formulate the field equations. Next, we shall solve these field equations using 
different ansatz of our parameters. Then we present the graphical representation of our 
results. Finally we conclude with the discussion of our results. 
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2 The model 



The static spherically symmetric spacetime (in geometrical units G = 1 
onwards) is taken as 



e-(-)dt' - e^^'Ur^ - r\d9^ + sin^ 9d(j)^). 
The Einstein field equations (EFE) for the above metric are 
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Here p is the energy density while Pr and pt are the radial and transverse pressure densities 
of the fiuid. The conformal killing equation (1) becomes 



L^gik — ii;k + ik;i — '^'Qik- 

The above equations give the following equations as 
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Integration of Eqs. (7-10) yield 
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where Cj, i = 1, 2, 3 are constants of integration. Making use of Eqs. (11-13) in (3-5), we 
can write 
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Figure 1: Plot for the variation of vs r( km 



Thus Eqs. (14-16) represent the EFE in terms of the conformal factor ip [9, 10, 14]. It 
is generally assumed that stars are the spherically symmetric anisotropic fluid, so one has 
to use the equations (2) - (6) for the initial consideration. Since we consider anisotropic 
star model admitting conformal motion , so all the parameters e^, e^, p, Pr, Pt could be 
found in terms of conformal factor ip . In other words, unless one knows the exact form 
of ip, one could not say any things ( i.e. all the physical properties of the anisotropic 
fluid ). Our model is new and interesting in the sense that we are the first authors who 
consider anisotropic stars with different form of anisotropy admitting conformal motion. 
The assumption of different anisotropy leads different differential equations of conformal 
factor ip i.e ip is constrained by the role of anisotropy. Recent observations on highly 
compact objects hke X ray pulsar Her X - 1 , etc support the existence of anisotropic stars. 
But it is still unknown the stars maintain what type of anisotropy. Since, conformal killing 
vector provides a deeper in sight in to the space time geometry, so we propose anisotropic 
star models with different anisotropy admitting conformal motion. In the next section, 
we shall solve these conformal field equations using different ansatz of the surface tension 
Pt — Pr- The surface tension actually determines the anisotropy in the stellar model. 



3 Solution of Conformal EFE 



Now we consider three different cases to get exact analytical solutions. Since the metric 
coefficient Qu is independent of conformal factor ip, so all the cases Qu assumes the same 
form. 



4 



3.1 Pt = npri < n < 1. 



The first case deals when pr and pt are Unearly proportional to each other. It is assumed 
that Pt is smaller then pr. By using above anisotropic relation, Eqs. (15) and (16) together 
will give a non-homogeneous and non-linear differential equation as 
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Solving this equation, we get 
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where B is the constant of integration and will be determined later. Hence we get the 
exact analytical form all the parameters as 
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Differentiation of Eq. (20) with respect to r yields the pressure gradient 
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Here, one can note that the pressure gradient is a decreasing function of r ( see figure 9). 
In order to get some physically meaningful solution, we need to determine the constant 
B. Since vanishing of radial pressure at the boundary is a consequence of the junction 
condition [15], so one can use 

p^(r = i?) = 0. (25) 

to find the constant B. 



Notice that Eq. (25) is motivated due to the fact that pressure gradient is a decreasing 
function of r. Thus using Eq. (20) in (25) yields 
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Figure 2: Plot for the variation of vs r( km ). Here we assume the radius of the star 
i? = 10 km and C3 = 1. 



Using Eq. (26) in (18), we get 
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The mass function is 
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Thus using Eq. (26) in (28) we obtain 
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Since for n > |, the solutions are inconsistent, so we neglect the cases n > |. We also 
observe that r = gives a singularity. However, the solutions are valid for some radius 
r > 0. Now, we calculate the subluminal sound speed 
as 
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The above expression should be less than one depending on the parameters. Thus sound 
speed does not exceed that of light as fulfillment of causality condition. 
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Figure 3: Plot for the variation of radial pressure Pr vs t{ km ). Here we assume the 
radius of the star R = 10 km and C3 = 1. 
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Figure 4: Plot for the variation of transverse pressure pt vs r( km ). Here we assume the 
radius of the star R = 10 km and C3 = 1. 
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Figure 5: Plot for the variation of energy density p vs r( km ). Here we assume the radius 
of the star = 10 km and C3 = 1. 




Figure 6: Plot for the variation of the conformal factor ip vs r( km ). Here we assume the 
radius of the star i? = 10 km and C3 = 1. 
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Figure 7: Plot for the variation of mass m(r) vs r( km ). Here we assume the radius of 
the star R = 10 km and C3 = 1. 




Figure 8: Plot for the variation of the pressure anisotropy pt — Pr vs r( km ). Here we 
assume the radius of the star R = 10 km and C3 = 1. 
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Figure 9: Plot for the variation of the gradient of radial pressure ^ vs r( km ). Here we 
assume the radius of the star R = 10 km and C'^ = 1. 



3.2 Pt — Pr = 8^(% + C2)? where ci and C2 are arbitrary constants. 



By using above anisotropic relation [7], we get from Eqs. (15) and (16) 

2ipip' 2^'^ (ci - 1) 
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Solving this equation, we get 
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where A = C|(ci — 1), B = C|c2 and Ki is an integration constant. Thus we 
exact analytical form of all the parameters as 
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Figure 10: Plot for the variation of e'^ vs r( km ). Here we assume the radius of the star 
i? = 10 km with C2 = .001 and C3 = 1. 




The pressure gradient is given by 
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Applying the similar procedure as done in the previous section, pr{r = R) =0 gives 
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As before, we calculate the subluminal sound speed ? I 'i's 1 = 1 ^ I 
as 

6(5r2lnr + fs:ir2-|)-2C|-3r2(251nr + 5 + 2i^i) 
' '~ r2(651nr + 55 + QKi) + {A - 6Br^ Inr - 2Br^ - QKir^ + 2C|) 

The above expression should be less than one depending on the parameters. Thus sound 
speed does not exceed that of light as fulfillment of causality condition. 
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Figure 11: Plot for the variation of radial pressure Pr vs r( km 
radius of the star R = 10 km with C2 = .001 and C3 = 1. 



Here we assume the 




Figure 12: Plot for the variation of transverse pressure pt vs r( km ). Here we assume the 
radius of the star i? = 10 km with C2 = .001 and C3 = 1. 
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Figure 13: Plot for the variation of energy density p vs r( km 
radius of the star = 10 km with C2 = .001 and C3 = 1. 
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Cj = -.l •Cj = -0.01— -Cj = -0.001 



Figure 15: Plot for the variation of mass m(r) vs r( km ). Here we assume the radius of 
the star R = 10 km with C2 = .001 and C3 = 1. 
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Figure 16: Plot for the variation of the gradient of radial pressure vs r( km ). Here 
we assume the radius of the star R = 10 km with C2 = .001 and C3 = 1. 
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Figure 17: Plot for the variation of the pressure anisotropy pt — Pr vs r( km ). Here we 
assume the radius of the star R = 10 km with C2 = .001 and C3 = 1. 

4 Conclusion and Discussion 

In the previous section, we have given the pictorial representation of the parameters in- 
volved in the three cases. The constant of integration is fixed at C3 = 1 while the radius 
of the compact star is chosen to he R = 10km. All the other physical parameters like 
the metric functions and e'^, radial pressure pr, transverse pressure pt, surface tension 
Pr — Pt, niass function m{r), density p(r) and the radial pressure gradient dp^/dr are 
plotted against the radial parameter r in all the cases. 

In figures 1-9, we have taken different values of n in the range (0,1) to determine param- 
eters from case-1. We observe that the radial and transverse pressures gradually decrease 
as r approaches to the surface R. This result is consistent with our boundary condition of 
vanishing pressures all at the stellar surface. It comes as no surprise that the pressure and 
the mass profile of the compact stars is all along similar to the normal hydrogen burning 
stars. More specifically, the mass profile of the star shows sharp increase for n = ^ 
while it is small for values n < j^. The surface tension pt — Pr of the star also decreases 
along with increasing r. Note that the surface tension varies along negative range due to 
higher radial pressure. Thus transverse pressure although exists but is smaller then the 
radial pressure along the allowed range of parameter r. The pressure gradient VrPr or 
dpr/dr also decreases outward along r. For small values of parameter n, the gradient is 
steeper then for large values. The density profile of the compact star falls for increasing 
r. The conformal factor ?/'(r) assumes maximum value at the star's center while it ap- 
proaches zero near the stellar surface. Therefore, the vanishing conformal factor at the 
star's surface could be used as an alternative boundary condition for solving conformal 
field equations. 

Similarly, the parameters obtained in case-2 are plotted in figures 10-17 using different 
values of parameter ci. The constant parameter is fixed at C2 = 0.001. The pressure 
profile of the star in this case is similar to the earlier in case-1 having steep slopes with 
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increasing r. Also, the density profile and the conformal factor show convergence for large 
values of r. The mass profile is steady for ci = —0.001 while it is steepest for ci = —0.9. 
Similarly the declining pressure gradient profile is observed for large r which is consistent 
with our boundary condition Eq.(25). 

In figures 18-25, we have made a comparison of various stellar parameters obtained in 
cases 1 and 2. It is observed that the metric function e^, conformal factor ip, radial pres- 
sure Pr and the transverse pressure pt of both cases converge as r tends to R. The mass 
parameter m(r) is steeper in the second case then the first one. Therefore the mass will 
increase fast radially in the linear model of the surface tension. Moreover, the surface 
tension profile of both models goes asymptotically parallel to each other. The density 
profile in the first case goes to zero while it remains constant in the second case in the 
asymptotic limit of large r. 

One can note that any static spherically symmetric space time admitting conformal mo- 
tion suffers r = singularity. In our model also, we had to tolerate r = singularity. 
Therefore the model does not exist at r = 0. However, our solutions are valid for some 
radius r > 0. We also note that there is a singularity in the mass density, in spite of, the 
total mass is finite. So, the models ( case - I and case - II ) are physically acceptable. In 
addition to this, one can find that the matter density and fluid pressure are non negative 
and gradient ^ is decreasing with r. Thus our models ( i.e. interior solutions of the 
gravitational field equations ) are fully physically meaningful. As our models suffer r = 
singularity, the minimum value of r is taken to be a positive quantity ( km ). Since radius 
of the star is taken to 10 km ( this is justified as the radius of some of the observed strange 
stars are nearly equal to 10 km ), so one can note that the total mass of the star is m ( 
r=10). Our results have been shown in fig. 7 and fig. 15. 

In our current analysis, we have used an idealization of spherical symmetry, stationary 
and static compact star. In general, from the astrophysical point of view, most of the 
observed or predicted compact star candidates are no longer static but rotating about a 
unique axis that may be oblique from the magnetic axis of the star. More exotic stars 
including magnetars (driven by mostly magnetic fields) and pulsars (mostly driven by 
higher angular momentum) are well known examples of rotating compact stars. In our 
forth coming work, we plan to work out a similar analysis presented here for the fast 
rotating or ultra-fast rotating compact stars. We also plan to work on the rotating stars 
with the slow rotation approximation as well. 
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Figure 18: Plot for the variation of vs r( km ). Solid line and dotted line represent 
the case 1 and case 2 respectively for suitable choices of the parameters. 
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Figure 19: Plot for the variation of conformal factor iIj vs r( km ). Solid line and dotted 
line represent the case 1 and case 2 respectively for suitable choices of the parameters. 



18 




123456789 10 



r 

Figure 20: Plot for the variation of energy density p vs r( km ). Solid line and dotted 
line represent the case 1 and case 2 respectively for suitable choices of the parameters. 
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Figure 21: Plot for the variation of mass m{r) vs r( km ). Solid line and dotted line 
represent the case 1 and case 2 respectively for suitable choices of the parameters. 
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Figure 22: Plot for the variation of the radial pressure pr vs r( km ). Solid line and dotted 
line represent the case 1 and case 2 respectively for suitable choices of the parameters. 
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Figure 23: Plot for the variation of the transverse pressure pt vs r(km). Solid line 
and dotted line represent the case 1 and case 2 respectively for suitable choices of the 
parameters. 
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Figure 24: Plot for the variation of the gradient of radial pressure vs r(km). Solid 
line and dotted line represent the case 1 and case 2 respectively for suitable choices of the 
parameters. 
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Figure 25: Plot for the variation of the pressure anisotropy pt — Pr vs t{ km ). Solid line 
and dotted line represent the case 1 and case 2 respectively for suitable choices of the 
parameters. 
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